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SECTION  1 


INTRODUCTION 


This  report  is  a  summary  of  work  performed  under  Con¬ 
tract  DNA001-79-C-0079 ,  entitled  "Precharge-Enhanced  SGEMP". 

The  aim  of  this  effort  was  to  develop  phenomenology  to  model 
the  enhancement  of  the  electromagnetic  response  of  spacecraft 
resulting  from  electrostatic  charging  in  the  natural  space  en¬ 
vironment  or  in  laboratory  simulations. 

It  is  now  well-accepted^  that  geosynchronous  spacecraft 
not  infrequently  experience  "magnetospheric  substorms"  which 
cause  them  to  attain  negative  potentials  of  'v-10  kv.  Further, 
spacecraft  can  charge  differentially  to  potentials  sufficient 
to  causa  spontaneous  electrostatic  discharges.  Commonly, 
these  discharges  cause  logic  malfunctions  requiring  correction 
by  ground-based  command.  In  a  few  instances  they  have  caused 
serious  damage  to  spacecraft. 

Another  well-studied  area  is  the  electromagnetic 

response  of  a  spacecraft  to  an  x-ray  burst  (SGEMP).  Typically, 

2 

an  x-ray  burst  causes  emission  of  ^100  A/m  of  electrons  with 
energy  of  'vl  keVr  together  with  larger  numbers  of  low  energy 

rol  _0 

electrons1  J  for  a  time  'vlO  seconds.  The  electromagnetic 
response  consists  of  surface  currents  'vl  A/m.  Normally,  the 
low  energy  electrons  play  a  negligible  role  in  this  response. 
Also,  increasing  the  x-ray  fluence  (and  thus  the  emitted  cur¬ 
rent)  does  not  greatly  increase  the  response  because  of  space 
charge  limiting  effects.  These  space  charge  effects  can  be 
somewhat  lessened  if  the  spacecraft  is  charged  negatively; 
however,  a  factor  of  two  is  the  largest  enhancement  expected 
from  such  uniform  charging. 

Much  larger  effects  are  expected  if  the  spacecraft  is 
differentially  charged.  The  basis  for  this  expectation  lies 
in  the  much  larger  energies  stored  in  differentially  charged 
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spacecraft.  We  at  Systems,  Science  and  Software ^  predicted 
SGEMP  enhancements  of  an  order  of  magnitude  under  these  condi¬ 
tions.  Such  effects  were  demonstrated  dramatically  in  measure- 
r  4 1 

ments  performed  by  Mission  Research  Corporation  (MRC)  on  the 
Skynet  Qualification  Model.  It  was  found  that  differentially 
precharging  the  model  by  electron  spraying  caused  SGEMP  en¬ 
hancements  of  a  factor  of  5.  Also  observed  was  DGEMP  — 
electromagnetic  response  due  to  spontaneous  discharges.  In 
one  ir reproducible  but  well-documented  case,  a  discharge  was 
triggered  by  an  SGEMP  pulse. 

These  findings  resulted  in  two  efforts  to  develop 
phenomenology  for  precharge-enhanced  SGEMP.  MRC  was  directed 
to  perform  experiments  using  simple  geometry  and  well- 
characterized  materials.  Systems,  Science  and  Software  was 
directed  to  perform  the  effort  described  in  this  report. 

Our  goals  were: 

1.  Provide  analytical  support  (using  the  NASCAP 
code) 15]  for  the  precharging  phase  of  the  MRC 
experiments. 

2.  Provide  consulting  services  on  the  general  ques¬ 
tion  of  laboratory  precharging  of  spacecraft, 
using  the  NASCAP  code  where  appropriate. 

3.  Perform  analysis  of  the  MRC  experiments. 

4.  Develop  analytical  tools  and  phenomenology  ap¬ 
propriate  to  precharge-enhanced  SGEMP,  electro¬ 
magnetic  response  to  discharges  and  related 
matters. 

In  Section  2  we  discuss  the  MRC  experiments  and  our 
support  of  them.  We  address  as  well  (Appendix  C)  the  question 
of  precharging  a  real  spacecraft  in  a  large  tank  by  electron 
spraying.  This  question  arose  in  connection  with  the  proposed 
Satellite  X-ray  Test  Facility  (SXTF) . ^  Analysis  of  precharge- 
enhanced  SGEMP  under  conditions  similar  to  the  MRC  experiments 
appears  in  Section  4  and  Appendix  B.  We  concluded  (as  did 
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MRCl'J,  in  accordance  with  our  previous  workl~'J)  that  the  en¬ 
hancement  was  associated  with  acceleration  of  low  energy 
electrons  by  strong,  localized  electric  fields. 


We  abstracted  two  generic  precharge  configurations  for 
theoretical  treatment.  The  first  (Appendix  A)  was  an  emitting 
surface  with  a  tangential  field.  This  led  to  the  concept  of 
the  conductivity  of  a  space  charge  layer,  as  well  as  a  pos¬ 
sible  mechanism  for  SGEMP-triggered  surface  flashover.  The 
conductivity  concept  later  proved  useful  to  estimate  the  skin 
current  in  the  metallization  beneath  the  emitting  dielectric 
(Section  4.3). 


The  second  generic  configuration  was  a  differentially 
charged  edge  (step  function  potential) .  This  configuration 
is  analyzed  in  Section  3.  A  first  surprising  result  (747 
effect)  was  that  the  anode  current  density  was  nearly  three 
orders  of  magnitude  below  the  cathode  current.  The  problem 
could  also  be  done  self-consistently  with  space  charge.  The 
"747  effect"  then  became  a  less  exciting  "38  effect",  and 
the  space  charge  limited  current  was  found  to  be  divergent. 
Also  in  Section  3  is  the  analytic  solution  for  the  potential 
of  a  charged  disk  on  a  ground  plane. 

Section  4  describes  a  new  hybrid  electrostatic/elec¬ 
trodynamic  2-D  code  designed  to  study  SGEMP  under  conditions 
where  high  local  resolution  is  required.  This  problem  has 
been  applied  both  ^o  the  MRC  type  problem  and  to  discharges, 
and  has  proven  to  be  an  excellent  analytical  tool. 

Finally,  Section  5  describes  a  theory  of  dielectric 
surface  discharges.  We  postulate  that  during  a  discharge 
the  entire  dielectric  surface  behaves  as  a  space  charge 
limited  emitter.  As  noted  above,  the  current  for  a  step- 
function  potential  is  divergent,  so  it  is  necessary  to  pos¬ 
tulate  a  voltage  profile  near  the  edge  of  the  dielectric. 

The  scaling  of  the  current  and  the  character  of  the  emission 


profile  is  treated  analytically  in  detail.  Additional  in¬ 
formation  comes  from  simulations  using  the  code  described  in 
Section  4.  We  find  that  this  theory  gives  good  agreement 
with  experimental  measurements. 

In  summary,  we  have  found  that  the  electromagnetic 
response  to  x-rays  under  precharged  conditions,  or  to  dis¬ 
charges,  is  dominated  by  the  flow  of  low  energy  electrons  in 
strong  local  fields.  We  have  proposed  a  mechanism  for  SGEMP- 
triggered  flashover  involving  these  low  energy  electrons.  We 
have  shown  how  amplification  of  the  SGEMP  response  by  an  order 
of  magnitude  over  the  non-precharged  response  is  possible.  We 
have  developed  both  computational  tools  and  approximate  theo¬ 
retical  models  to  aid  in  predicting  this  response.  Such  models 
could  be  incorporated  as  "effective  source"  treatments  in  a 
practical  3-D  SGEMP  code.  In  addition,  we  have  shown  NASCAP 
to  be  a  helpful  tool  for  analysis  of  precharging  under  labora¬ 
tory  conditions  and  for  design  of  experiments. 


SECTION  2 


EXPERIMENTAL  SIMULATION*  OF  PRECHARGE- ENHANCED  SGEHP 

The  experimental  program  designee  to  study  precharge- 
enhanced  SGEHP  is  described  in  this  section.  A  very  brief 
outline  of  the  experimental  details  is  given  in  Section  2.1. 
NASCAP  charging  simulations  performed  by  S^  to  aid  in  the 
experimental  design  and  setup  are  presented  in  Section  2.2. 

The  main  results  of  the  experiments  are  summarized  and 
assessed  in  Section  2.3. 

2.1  DESCRIPTION  OF  EXPERIMENTS 

During  the  early  part  of  1979  experiments  were  per¬ 
formed  by  Mission  Research  Corporation  (HRC)  using  the  OWL  II* 
photon  source  to  study  precharge-enhanced  SGEHP.  KRC  has  de¬ 
scribed  these  experiments  elsewhere. ^  This  discussion  is 
included  for  completeness,  and  to  provide  a  context  for  the 
analytical  and  simulation  work  which  follows. 

The  basic  experiment  consisted  of  an  aluminum  cylinder, 

1.2  m  in  diameter  and  51  cm  long,  with  an  £2  cm  diameter  hole 
in  one  end  against  which  a  sample  was  placed-  we  will  be 
concerned  with  two  of  the  four  samples  used:  0.013  cm 
silvered  teflon  mounted  on  a  copper  backplate,  and  0-05  cm 
white  paint  on  fiberglas  mounted  on  a  segmented  backplate. 

The  experiment  was  located  in  an  approximately  cylindrical 
tank  with  effective  radius  "-1.2  m  and  length  "-3.5  ~. 

The  insulating  surfaces  were  charged  by  an  electron  gun 
-9  2 

with  ”vl0  A/cm  current  over  a  period  of  several  minutes  to 
potentials  of  10  kV  to  15  k V .  The  potential  was  far  from 
uniform  along  the  surface;  it  was  commonly  measured  to  vary 
several  kilovolts  on  a  few  centimeter  scale. 
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The  primary  measurement  taken  upon  photon  exposure 
was  the  current  between  the  aluminum  cylinder  and  the  sample 
backplate  (teflon  sample) ,  or  the  currents  between  backplate 
segments  (paint  sample) .  Photodiodes  were  used  to  determine 
pulse  characteristics,  and  the  final  object  potential  was 
measured.  During  the  time  of  the  photon  pulse  the  cylinder 
was  effectively  ungrounded. 

2.2  NASCAP  PRECHARGING  CALCULATIONS 

The  NASCAP  code  was  used  to  simulate  the  precharging 
of  the  MRC  test  object.  The  results  of  these  calculations 
were  described  to  MRC  experimenters  in  advance  of  their  SGEMP 
test  series  to  assist  them  in  the  design  and  interpretation 
of  their  experiments. 

NASCAP  allows  a  detailed  simulation  of  the  charging 

of  an  obj  .ct  in  a  test  tank  environment,  as  documented  else- 
f  si 

where.  A  charging  surface  is  described  in  terms  of  its 

thickness,  dielectric  constant,  energy  and  angle  dependent 
emission,  as  well  as  bulk  and  surface  conductivity.  The 
initial  beam  current  profile  can  be  specified  at  the  sample 
surface,  and  electrostatic  effects  on  the  incident  electron 
trajectories  are  included  as  the  charging  proceeds.  Compli¬ 
cated  current  profiles,  such  as  that  produced  by  a  rastered 
beam,  can  be  included  in  the  simulation. 

The  NASCAP  model  of  the  SGEMP  test  object  is  illus¬ 
trated  in  Figures  2.1  and  2.2.  NASCAP  solves  Poisson's 
equation  on  a  cubic  mesh,  on  which  an  octagonal  prism  was 
used  to  model  the  aluminum  cylinder.  The  insulating  surface 
was  normal  to  the  axis  of  the  electron  gun,  located  approxi¬ 
mately  2  m  above  the  test  object.  The  cylinder  and  the 
surrounding  tank  were  grounded  throughout  the  charging  simu¬ 
lations. 
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Figure  2.1.  NASCAP  model  of  cylindrical  SGEMP  test  object 
Zone  size  =  8.6  cm. 


MATERIAL  lEoEND 


Figure 


Results  ,of  the  charging  study  for  a  teflon  insulating 
surface  using  a  10  kV  beam  are  illustrated  in  Figures  2.3 
through  2.5.  The  initial  electron  beam  profile  was  taken 
to  be  uniform  across  a  circular  region  of  variable  radius, 
while  the  sample  radius  was  fixed  at  40  cm.  The  MRC  experi¬ 
menters  expressed  an  interest  in  predicted  voltage  profiles 
as  the  initial  radius  of  the  charging  beam  was  reduced  below 
the  insulating  sample  radius,  so  that  the  edges  of  the  sample 
might  remain  near  zero  potential.  The  voltage  profile  for 
an  initial  beam  radius  of  40  cm  was  almost  indistinguishable 
from  that  for  the  30  cm  case,  illustrated  in  Figure  2.3.  The 
strong  electric  fields  near  the  edge  of  the  sample  deflect 
the  30  cm  beam  so  that  the  entire  surface  is  sprayed  with 
electrons.  The  potential  contours  during  charging  which 
cause  the  beam  spreading  are  shown  in  Figure  2.4,  and  typical 
trajectories  are  shown  in  Figure  2.5a.  In  both  cases,  the 
sample  charges  less  near  the  edges  due  to  a  higher  yield  of 
secondary  electrons  for  non-normal  incidence.  Also  shown  in 
Figure  2.3  is  the  voltage  profile  resulting  from  a  beam  ini¬ 
tially  10  cm  in  radius.  In  this  case  the  beam  spreading, 
shown  in  Figure  2.5b,  is  not  as  large  since  the  integrated 
deflecting  forces  are  much  reduced  for  such  a  collimated  beam. 
A  significant  fraction  of  the  insulator  remains  uncharged, 
and  the  potential  gradient  is  large  near  the  beam  edge.  As 
a  result  of  the  above  study,  MRC  experimenters  were  aware 
before  their  test  series  that  sample  shields  or  very  re¬ 
stricted  beams  would  be  required  if  the  effects  of  a  sharp 
potential  gradient  located  some  distance  from  the  sample 
edge  were  to  be  investigated. 
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SURFACE  POTENTIAL 
(VOLTS) 


Figure  2.3.  Voltage  profile 
voltage  =  10  k\ 


iMPLE  CENTER  (CM) 
across  teflon  surface 


during  charging  of  teflon 


Figure  2.5a.  Particle  trajectories  at  equilibrium.  Initial 
beam  radius  =  30  cm. 


Figure  2.5b.  Particle  trajectories  at  equilibrium.  Initial 
beam  radius  =  10  cm. 


2.3  EXPERIMENTAL  RESULTS  AND  CRITIQUE 

For  comparison  with  results  to  be  presented  later, 
there  were  three  main  results  from  these  experiments: 

1.  Increased  Surface  Currents 

Under  precharged  conditions  surface  currents  of 
10-20  A/m  were  measured  at  the  sample  edges.  By 
comparison,  non-precharge  shots  had  surface  cur¬ 
rents  an  order  of  magnitude  less. 

2.  Lack  of  Surface  Discharging 

In  general,  no  significant  dielectric  discharging 
was  seen.  This  was  expected:  given  10  pf/cm2  and 
10  mA/cm2  for  10“7  seconds,  AV/V  =  100  volts/ 

15  kilovolts.  However,  in  one  or  two  cases  sub¬ 
stantial  local  redistribution  of  charge  was  ob¬ 
served. 

3.  Blowoff  Charge 

The  potential  reached  by  the  body  is  a  measure  of 
the  blowoff  charge.  Typically  this  potential 
was  1-2  kV,  indicating  a  blowoff  of  ^0.3  yC.  In 
at  least  one  case,  however,  an  additional  ^1  yC 
of  charge  was  blown  off  during  ,vl  ysec  following 
the  shot.  MRC  referred  to  this  phenomenon  as 
"delayed  photocurxent" . 

Most  of  the  difficulties  with  these  experiments  were 
beyond  the  control  of  MRC.  The  small  number  of  photon  shots 
and  their  irreproducibility  make  the  data  difficult  to  inter¬ 
pret.  There  were  questions  of  debris  being  shot  into  the 
experiment  area.  Also,  the  nonuniform  potentials  on  the 
dielectric  surfaces  complicates  some  interpretations.  Be¬ 
cause  of  these  factors,  planned  experiments  such  as  narrow 
or  patterned  charging  profiles  had  to  be  abandoned. 

As  theorists  attempting  to  obtain  data  for  interpre¬ 
tation  we  found  that  the  data  taken  as  Polaroid  scope  pictures 
was  very  difficult  to  work  with.  These  data  included  the 
object  potential  measurement  and  the  back-biased  diode  data 
used  to  characterize  the  photon  pulse.  The  diodes  having 
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insulating  materials  were  in  unknown  states  of  precharge  and 
gave  conflicting  results  from  shot  to  shot.  Also  these  pictures 
reproduced  so  badly  as  to  be  often  unreadable  when  we  obtained 
them. 

Despite  the  above  difficulties,  we  believe  these  ex¬ 
periments  provide  good  estimates  of  the  amount  of  SGEMP  en¬ 
hancement  expected  from  precharging.  They  can  be  used  to 
validate  predictive  theories  and  computer  codes.  The  delayed 
effects  which  were  seen  represent  a  type  of  discharge,  and 
show  that  nuclear  photons  can  trigger  weak,  localized  dis¬ 
charges.  Even  though  strong  triggered  discharges,  such  as 
observed  in  the  previous  Skynet  tests,  were  not  seen  here,  we 
feel  that  such  strong  synergistic  effects  remain  a  real  pos¬ 
sibility. 
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SECTION  3 


ANALYTICAL  CALCULATIONS 


This  section  describes  the  analytical  solution  of 
three  model  problems  which  illustrate  the  behavior  of  electro¬ 
static  potentials  and  electron  trajectories  near  sudden  sur¬ 
face  potential  changes.  Although  these  model  problems  are 
considerably  cinplified  abstractions  of  real  experiments, 
they  provide  insight  into  the  central  physical  mechanisms, 
and  they  generate  parametic  representations  of  SGEMP  response 
which  can  represent  or  even  replace  the  results  of  more  de¬ 
tailed  simulations.  Potentials  and  electron  trajectories 
near  the  interface  of  infinite  half  planes  is  discussed  in 
Section  3.1,  and  the  space  charge  limited  flow  of  electrons 
between  such  half  planes  is  analyzed  in  Section  3.2.  The 
analytic  potentials  around  a  charged  disc  are  presented  in 
Section  3.3. 

3.1  ANALYTIC  TREATMENT  OF  ELECTROSTATIC  POTENTIAL  AND 

PARTICLE  TRAJECTORIES  NEAR  SUDDEN  SURFACE  POTENTIAL 
CHANGES 

The  largest  effects  of  differential  charging  in  SGEMP 
occur  in  the  immediate  vicinity  of  edges  between  charged  and 
uncharged  surfaces  (e.g.,  charged  dielectric  and  exposed, 
grounded  conductor) .  Naively,  one  would  expect  large  cur¬ 
rents  of  electrons  to  follow  field  lines  from  one  surface 
to  the  other.  This  turns  out  not  to  be  the  case;  in  fact, 
the  current  density  incident  on  the  collecting  surface  is 
nearly  three  orders  of  magnitude  less  than  the  current  density 
on  the  emitting  surface. 

Consider  the  situation  shown  in  Figure  3.1.  Let  r 
denote  the  radial  distance  from  the  edge,  and  6  the  angle, 

(see  Figure  3.2).  The  electrostatic  potential  <j>(r,9)  satis¬ 
fies  the  Laplace  equation 


Figure  3.2.  Geometrical  construction  for  particle  trajectories  and 
space  charge  density  near  charged  edge. 


v2$ (r,0) 


0 


(3.1) 


with  boundary  conditions 


4>(r ,0=u)  =  +Vq/  <?  (r ,0=0)  =  0 


(3.2) 


where  r  is  the  distance  from  the  edge  and  8  the  angle  measured 
from  the  grounded  to  the  positive  surface.  The  potential 
<?(r,0)  which  satisfies  (3.1)  and  (3.2)  is: 


<P  (r,d)  =  +VQ  9 /if 


(3.3) 


Electrons  emitted  from  the  positive  surface  will  be  immediately 
limited  by  the  electrostatic  field,  so  we  need  to  consider  emis¬ 
sion  only  from  the  grounded  surface. 

Let  us  find  the  trajectory  of  an  electron  leaving  from 
the  left  half  plane,  with  initial  position  (rQ,»).  The 
Lagrangian  is  given  by: 

•2' ( r , 8 )  =  T-V  =  T+e4>  =  4  n»r2  +  i  mr2'*2  +  eV  2-  (3.4) 

i  z  o  *> 


It  is  advantageous  to  write  &  in  terms  of  the  dimensionless 
radius,  R  =  r/r Q,  and  time 


T  _  t_ 
r 


(!!o) 

\mTi  / 


1/2 


(3.5) 


In  terms  of  R,  t  and  using  primes  to  denote  the  derivative 
with  respect  to  t  one  obtains 


irrj,=  i  (R'2  +  r29,2)  ' 

o 


(3.6) 


Thus,  if  initially  (x  =  0)  the  dimensionless  angular  momentum 
2 

R  9 ’  is  zero,  it  follows  that 
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(3.7) 


and 


R3R* ’  =  t2  (3.8) 

Equations  (3.7)  and  (3.8)  together  with  energy  conser¬ 
vation,  result  in  a  single  dimensionless  particle  trajectory 
for  electrons  emitted  with  negligible  initial  energy 
{R‘ (x=0)  =  3’ (t=0)  =  0] : 

12  2  2 

tJ  =  +  |  (R*Z  +  Tz/R  )  (3.9) 

We  have  solved  for  the  single  trajectory  numerically;  we  show 
it  in  Figure  3.3.  Defining  (X,Y)  =  (Rcos3,  RsinS) ,  this  tra¬ 
jectory  has  the  following  properties: 

•  The  particle  returns  back  to  the  plane  on  the  left 
hand  side,  (3=2),  of  the  edge  at  X  =  R  i  -747,  at 
x  =  314. 

•  Y  =48  and  occurs  at  X  =  250. 

max 

•  Y(x=0)  *  2.5. 

•  dY/dX  (6=2)  -0.17  (angle  of  incidence  =  80°). 

It  is  instructive,  as  an  example,  to  notice  that  at  VQ  =  10  kV 

and  r  =  1  cm  the  time  taken  to  return  to  the  plane  is  132  * 

-9  ° 

10  sec.  We  may  conclude  that,  in  the  absence  of  space 
charge  effects,  current  emitted  from  a  negatively  charged 
surface  returns  to  a  neighboring  position  only  after  a  long 
time  and  at  low  current  density  (down  by  747).  Effectively, 
this  current  may  be  considered  to  escape  the  neighborhood  of 
the  spacecraft  or  test  object. 

We  now  can  derive  the  charge  density  of  electrons 
moving  in  this  Laplace  field.  We  can  do  that  analytically 
because  of  the  existence  of  a  unique  trajectory.  Consider 
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Figure  3.3.  Particle  trajectory  near  charged  edge:  (a)  full  trajectory 
(b)  initial  portion  of  trajectory. 


electrons  emitted  at  r  <  r  <  r  +  dr  and  between  t  and 

o  o  o 

t  +  dt.  Conservation  of  charge  implies 


ffe1  dr  ae  -ffJo 


dt  dr 


(3.10 


(See  Figure  3.2.)  Since  Eq.  (3.10)  holds  for  an  arbitrary 
domain  of  integration  it  follows  that 


P  =  J 


3(6, r) 


3(t,r0) 


(3.11) 


where  3 (0,r)/3 (t,r  )  stands  for  the  Jacobian  of  (0,r)  with 


respect  to  (t ,r  ) . 


Using  the  scaling  implied  by  (3.5)  and 
Eqs.  (3.7)  and  (3.9)  we  obtain 


P(r,6)  = 


(fo) 

\  irm  / 


1/2 


(3.12) 


where  x  is  the  dimensionless  time  it  takes  an  electron  to 
reach  point  (r,0) . 

Equation  (3.12)  allows  us  to  obtain  (non-self-consist- 
ently)  the  space  charge  limiting  field  as  a  function  of  time 
and  distance  from  the  edge.  This  should  provide  a  simple  way 
of  estimating  tne  SGEMP  enhancement  due  to  differential 
charging. 
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3.2  ANALYTIC  TREATMENT  OF  SPACE  CHARGE-LIMITED  EMISSION 
NEAR  SUDDEN  SURFACE  POTENTIAL  CHANGES 

Problem  Definition,  Geometry  and  Boundary  Conditions 

Consider  the  half  space  Z  _>  0 ,  bounded  at  Z  =  0  by 
an  infinite  plane,  the  left  half  (0  =  tt)  is  at  potential 
V  ,  right  half  plane  (0  =  0)  is  grounded,  electrons  are 
emitted  with  zero  initial  velocity,  and  <|>(r,o)  0.  What 

is  the  space  charge-limited  emission  j(r,o)  from  9=0? 


Z 


Dimensional  Considerations 

The  solution  of  Poisson's  equation  and  Newton's  equa¬ 
tions  consistent  with  the  boundary  conditions,  and  the  initial 
conditions  of  electron  motion,  can  be  expressed  in  terms  of 

the  dimensional  quantities  r  ,  t,  V,  e,  and  m,  where  r  is 

o  o 

the  initial  radius.  We  consider  these  to  be  the  independent 
variables  of  the  problem.  The  only  dimensionless  combination 
that  can  be  formed  from  these  quantities  is 
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where  v  is  given  by 


1  2 

2  mvo  =  eV 


Consequently,  we  may  write: 


Electron  radial  and  angular  positions 


r  =  rQ  R(r) 


9  =  0(t) 


Electric  potential 


<P  =  V  <4>  (  t  ) 


Electron  density 


n  =  — j  n(t) 


Electron  particle  current  density 


1  /2e\1/2  V3//2  , 

3=2 ?(i-)  —  J(T) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


(3.18) 


(3.19) 


Observe  in  particular  that  the  solutions  are  separable  into 
products  of  functions  of  rQ  and  functions  of  x.  Because  of 
Eq.  (3.16),  this  separability  is  preserved  if  9  rather  than 
x  is  used  as  an  independent  variable. 


Solutions 


The  remainder  of  this  section  will  be  directed  toward 
obtaining  R(t)  ,  0(x),  N(x)  and  J(x).  Of  particular 

interest  is  the  current  density  on  the  emitting  surface. 

Consider  first  the  equations  of  motion 


m 


(3.20) 


(3.21) 


where 


With  the  help  of  Eas.  (3.13)  through  (3.17),  these  equations 
become 


R(R"  -  RO'2) 


(3.22) 


where 


subject  to  the  initial  conditions 


(3.23) 


0'  =  0  at  x  =  0  (3.24) 

Later,  in  analyzing  the  potential  equations  we  will  require 
transformations  from  r,9  to  rQ,T: 
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(3.25) 


From  Eqs .  (3.15)  and  (3.16) 


Now  consider  Poisson’s  equation 


I  L  r  M  +  L  32<f>  _ 
r  3r  r  3r  r2  ^2  " 


47ine 


(3.26) 

(3.27) 

(3.23) 

(3.29) 

(3.39) 

(3.31) 

(3.32) 

(3.33) 
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In  steady  state  the  electron  density  is  related  to  the  dynamical 
variables  by  the  conditions  of  charge  conservation 


/ 

V (r, 0) 


nr  dr  d6 


t) 


j(r„,0)  dr  dt 


(3.34) 


where  the  region  V  in  r,9  space  is  determined  from  the  emis¬ 
sion  region  S  in  r  , t  space  by  that  dynamical  transformation 
defined  by  the  equations  of  motion.  Since  we  may  choose  S 
arbitrarily 


or 


nr 


3  (r,  0) 


3(Vt) 


=  j(ro,0) 


nr 


f/ar  \  1 

'M\  -  ll 

)r\  /  36  \ 

_ 

0 

M 

CD 

,9t/r  v 

O 

5t/r  \3ro/, 
o  t 

=  3(ro,0) 


(3.35) 


(3.36) 


From  (3.15)  and  (3.16) 


whence 


j(ro'0) 

v  R20' 
o 


(3.41) 


Return  now  to  Poisson's  equation  and  evaluate  the  radial 
derivatives 


i_  r  M  =  JL  _i_  /r  \ 

3r  3r  roR  3rQ  [ro  l  8r 


(3.42) 


where  on  the  right  hand  side  3/3rQ  is  taken  at  constant  t 
and  on  the  left  hand  side  3/Dr  is  taken  at  constant  0. 

Using  Ea.  (3.22) 

?  h  r  H  =  FT  (£-  E(R"  -  r®'2>)  ■  0 

O  \  O  T 

since  R  and  G  are  functions  of  x  alone.  Similarly  we  find 

1  =  2V  32R2G"  , 

r  392  0'  dx2 


(3.43. 


(3.44) 


Recalling  Ecs.  (3.19)  and  (3.41),  Poisson's  equation 


reduces  to 


Isifi  .  J(0) 

dxz 


(3.45) 


From  Eq.  (3.21)  and  the  required  vanishing  of  initial  veloc¬ 
ity  and  normal  electric  field  at  6  =  0,  t:ia  solution  of  (3.45) 


(3.46) 
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Now  from  Eqs.  (3.22)  and  (3.31) 


R2(R"-R0'2)  =  /|i_\  .  K  (3.47) 

-  '"o/T 

where  K  is  a  constant  independent  of  r  and  x.  Since  d>  must 

o 

vanish  at  rQ  •*  °°  on  0  =  0  (x=0) 

K  =  0  (3.48) 

so  that 


$(rQ,0)  =  0 


(3.49) 


R3R~  =  i  J2  (0)  t4 


(3.50) 


The  constant  J(0)  follows  from  energy  conservation  and  the 

boundary  condition  $  =  1  on  the  reduced  potential  on  the  J 

surface  0  =  tt  where  x  =  x  : 

°  i 

R  2  (xq)  +  R2(to)20'2(xo)  =  1  (3.51)  j 

i 

i 

! 

The  reduced  emission  current  J(0)  is  the  value  that 
renders  the  solutions  of  (3.46)  and  (3.50)  consistent  with 
the  condition  (3.51).  The  evaluation  of  J(0)  is  simplified 
by  the  substitution 


s  =  J(0)1/3  x  =  J(0)1/3  — - 

o 


giving 


(3.52) 


,3^,1  4 

as2  < 


(3.53) 


2  dG  _  1  2 

ds  2 


(3.54) 
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=  J(0) 


(3.55 


with 


(aS  (so>)  -  R 

s  determined  from 


0 


ai  «•-> 


o») 


/ 


o  2 . 
s  ds 


(3.56 


The  solution  is  trivially  generalized  to  the  case  of  an  angle 
a  7^  n  between  the  two  half  planes. 

Particle  Trajectory 

The  problem  reduces  to  the  determination  of  the  unique 
particle  trajectory  established  by  the  functions  R  and  0  (Eqs. 
3.15,  3.16)  which  are  solutions  of  Eqs.  (3.53),  3.54). 

The  constant  J (0 )  is  given  by  Eq.  (3.55).  The  space 
charge  limited  electron  current  density  on  the  half  plane 
emitting  surface,  j  ,  is  related  to  J(0)  by 


j  =  1_  (  2e\1/2  J(0) 

3HP  2tt  \  m  /  2  J(0)  (3.57 

/  r 

o 

The  trajectory  was  generated  by  numerical  integration 
of  Eqs.  (3.53)  and  (3.54)  using  a  second-order  leapfrog  scheme, 
with  boundary  conditions  R(0)  =  1,  R'(0)  =  0,  0(0)  =  0.  Fig¬ 
ure  3.4  displays  the  solution.  An  electron  emitted  at  r  =  1 
from  the  9=0  half-plane  rises  to  a  maximum  height  of  6.7 
before  returning  to  the  9  =  »  half-plane  at  r  =  38.  (The 
corresponding  figures  for  the  zero  space  charge  case  are 
maximum  height  of  48  and  return  at  r  =  747.)  For  V  =  -10  KV 
and  r  =  1  cm,  the  time  to  return  to  the  plane  is  11  nanoseconds 


-«.•  -B.0  -30.0  -a.a  -30.0  -15.0  -10.0  -4.0  0.0  5.0 

X-flXIS 

Unique  particle  trajectory  from  solution  of  equc- 
tions.  Electrons  start  on  the  <p  =  0  half  plane 
(positive  X-axis)  at  X  =  1.0  and  return  to  the 
<j>  =  V  half  plane  (negative  X-axis)  at  X  =  -38.2. 
Note  disparity  in  the  X  and  Y  axis  scales. 


Emission  Current  Density 

For  half-planes  separated  by  an  angle  a 
J(0)  =  8.6  x  10  Therefore 


—  ii » 


jHP  -  9.0  x  10 


-8  V 


3/2 


(3.58) 


-2  .  _  . 

where  V  is  in  volts,  and  will  be  m  amps-cm  if  rQ  is  m 
cm.  The  parameter  J{0)  has  also  been  calculated  as  a  function 
of  a,  the  angle  between  half-planes,  and  the  results  are  dis¬ 
played  in  column  1  of  Table  3.1.  The  radial  distance  from  the 
origin  at  the  end  of  the  trajectory  is  given  in  column  2  of 
the  same  table. 

The  above  results  for  emission  current  density  are 
surprisingly  close  to  a  simple  estimate  obtained  from  the 
Child-Langmuir  expression  for  space  charge  limited  flow  be¬ 
tween  parallel  plates.  If  the  distance  between  plates  is  d, 
the  Child-Langmuir  current  density,  jCL,  is 


.  _  1_  (2e\ 

3CL  9-  \m  } 


1/2  v3/2 


(3.59) 


If  the  distance  d  is  estimated  by  d  =  arQ,  we  have 


JHP  9  T/_,  2 
rj —  =  j  J(0)g 

3cl 


(3.60) 


The  right  hand  side  of  Eq.  (3.60)  is  given  in  column  3  in 

Table  3.1.  As  expected,  the  ratio  approaches  unity  for  small 

angles.  At  large  angles,  the  ratio  deviates  surprisingly 

little  from  unity,  even  though  an  inspection  of  the  final  r 

values  in  column  2  indicates  that  the  estimate  d  =  ar  is 

o 

grossly  in  error.  However,  the  space  charge  limiting  is 
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certainly  much  less  in  the  two-dimensional  problem  than  in 
the  Child-Langmuir  case,  apparently  nearly  compensating  for 
the  error  in  the  estimate  of  d. 


Table  3.1.  Numerical  Results  for  Space  Charge  Limited  Flow 
Between  Half  Planes 


a 

J(0) 

■A. 

^hpA 

0.1 

22.2 

1.003 

1-00 

0.5 

0.88 

1.08 

0.99 

1.0 

A  0*1 

W  •  4.X 

1  T 

X.  J  / 

0.93 

m/2 

0.080 

2.29 

0.89 

2.0 

0.044 

4.16 

0-79 

2.5 

0.023 

10.5 

0.65 

3.0 

0.011 

28-9 

0.44 

7T 

0.0086 

38.2 

0.38 

a  =  angle  between  half  planes  in  radians,  J(0)  = 
unitless  parameter  determined  from  Eq.  (6),  r= 
distance  of  particle  from  origin,  jap/icL  =  rat-i° 
of  calculated  current  to  an  estimated  current 
using  the  Child-Langmuir  expression. 


3.3  ANALYTIC  SOLUTION  FOR  TIE:  ELECTROSTATIC  POTENTIAL  OF 
A  CHARGED  DISK  ON  A  GROUND  PLANE 

A  more  realistic,  though  less  analytically  tractable, 
case  is  that  of  the  electrified  disk  of  radius  rQ,  charged 
up  to  unit  voltage,  and  surrounded  by  a  grounded  plane - 


The  potential  upon  neglecting  space  charge  satisfies 
Laplace's  equation.  In  cylindrical  coordinates  p  i  r/r^. 
and  z  =  Z/r _,  we  have: 
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2  >  0 


(3.61) 


(P/Z) 


0 


}  (d>1,  z=0)  =  1  (3.62) 

6(p>l.  z=0)  =  0  (3.63) 

The  integral  representation  of  the  solution  for  (3.61- 
[81 

3.63}  is  known,  '  i.e.. 


,  %  2  T 

(p ,  z )  =  -  —  Iia 


a 

/ 


t  dt 


0  Vl-t' 


2  2 
>  +  (z+it) 


1/2 


(3-64) 


Introducing  the  substitution  t  -  sin©,  ©(p,z)  is  found  to  be 


~/2 

$(P,z'j  =  |  a© 


sinO  sin  | 

f i  tan  *  j 
[2  1 

f  2z  sin®  \ 

1 

U2  +  z2  -  sin-* 

!o/, 

,2  2 
(?  -*  z 

-  sin2©} 

2  2 
+  (2z  sinO) 

1/4 

(3.65) 


In  Figure  3-5  we  show  a  contour  plot  of  (;,z)  .  It 

is  worth  remarking  that  if  decays  faster  or  as  fast  as 
2  2  -1 

(p  +  z  )  and  that  the  monopole  and  dipole  moments  of  the 
induced  charged  density  are  zero.  No  higher  order  multi- 
poles  exist . 

Khile  this  case  is  far  more  complex.  It  would  be  use¬ 
ful  to  study  particle  trajectories  and  space  charge  for  both 
the  cases  of  negatively  charged  dielectric  on  a  ground  plane, 
and  for  a  positive  Interconnect-  The  delay  In  space  charge 
limiting  relative  to  the  limiting  time  for  an  uncharged  sur- 
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a  simple  estimate  of  SGEMP  enhancement. 


SECTION  4 


A  HYBRID  ELECTROSTATIC/ELECTRODYNAMIC 
SGEMP  COMPUTER  CODE 

4.1  WHY  DEVELOP  YET  ANOTHER  2-D  SGEMP  CODE? 

A  major  effort  under  this  contract  has  been  develop¬ 
ment  of  a  new  SGEMP  code  appropriate  to  spacecraft  on  test 
bodies  having  strong  differential  precharge.  The  special 
problem  here  is  that  electrostatic  potentials  vary  rapidly 

in  space  on  a  scale  of  1  cm  or  less.  In  fact,  one  would 

-2 

like  to  resolve  layers  of  dielectric  ^10  cm  thick.  A 
standard  SGEMP  code  with  uniform  1  cm  zoning  would  require 

4 

10  zones  per  square  meter  of  mesh  space,  and  a  timestep 

(At  <  Ax/c/2)  of  ^10  ^  seconds.  Yet  the  shortest  time  of 
-9 

interest  is  'vlO  seconds.  Nonuniform  zoning  with  otherwise 
standard  techniques  generates  spurious  electromagnetic  noise 
everywhere,  but  especially  at  the  irregular  interfaces.  A 
new  approach  was  required  to  achieve  a  practical,  generaliz- 
able  method  for  this  class  of  problem. 

Our  solution  was  to  solve  separately  the  conflicting 
parts  of  the  problem,  viz.,  high-resolution  electrostatics 
coupled  to  coarse-resolution  electrodynamics  by  the  tra¬ 
jectories  of  emitted  electrons.  Separation  of  fields  into 
electrostatic  (longitudinal,  curl-free)  and  electromagnetic 
(transverse,  divergence-free)  components  is  a  well-known 
technique.  For  the  problems  of  interest  the  longitudinal 
fields  are  the  primary  influence  on  the  particle  trajec¬ 
tories,  while  the  transverse  fields  determine  the  SGEMP 
response. 

In  using  this  code  to  model  the  MRC  test  body,  we 
resolved  the  f iberglas-paint  dielectric  layer  with  thickness 
0.02  cm,  had  spatial  zones  in  front  of  the  surface  with 
thickness  0.2  cm,  and  resolution  as  small  as  1  cm  along  the 
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dielectric  surface  near  the  edge.  (The  latter  could  have  been 
made  still  smaller  without  severely  testing  our  storage  re¬ 
strictions.)  The  electrostatic  grid  is  shown  in  Figure  4.1. 

The  electromagnetic  grid  (Figure  4.2)  had  uniform  6.67  cm 

-9 

zones,  and  a  timestep  of  0.125  x  10  seconds  was  used.  (Typ¬ 
ically  the  electrostatic-particle  portions  of  the  code  need 
be  invoked  only  once  every  2-8  timesteps.)  Excellent  results 
were  obtained. 


4 . 2  CODE  DESCRIPTION 


4.2.1  Mathematical  Summary 

The  electrostatic  portion  of  the  code  calculates  the 
electrostatic  potential,  V,  from 


7  v  =  -p/e. 


(4.1) 


and  tracks  particles  in  the  electric  field  E  =  E1  -  VV.  The 
electromagnetic  code  uses  the  particle  information  to  generate 
a  source  current  density,  J,  and  timesteps 


9r  /i  \ 

—  =  —  VXB-J  +  V^ 

o  9 1  \vo  -  ~  ~  -  v) 


(4.2) 


3B 

^  =  -V  X  EJ 
at 


(4.3) 


where  ij;  is  determined  by  requiring  V  •  E  =  0 ,  or 


-V  ty  =  V 


(  ~  v 


X  B  -  J 


(4.4) 


subject  to  the  same  boundary  conditions  as  V.  Knowledge  of 
B  is  then  equivalent  to  knowledge  of  surface  currents. 
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Figure  4.1.  Electrostatic  grid  used  for  precharge-enhanced 
SGEMP/discharge  calculations.  Note  fine  resolu¬ 
tion  near  dielectric-metal  edge  and  in  front  of 
emitting  surface. 
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The  centering  for  the  electromagnetic  code  is  shown  in 

Figure  4.3.  This  centering  has  the  advantages  that  currents 

and  potentials  are  defined  on  the  object  surfaces,  while 

electric  fields  are  defined  in  space,  making  it  easy  to  feed 

electromagnetic  fields  into  the  electrostatic  code  and  track 

particles  in  the  total  field  E  =  +  E^.  (For  problems  con- 

.  T  ~ 

sxdered  here,  E  is  a  small  correction.)  It  is  essential  that 
the  discrete  representations  of  Laplacian,  divergence,  and 
gradient  form  a  consistent  set,  so  that  V  •  E  exactly  van¬ 
ishes  numerically.  For  this  centering,  it  is  not  possible 
to  construct  a  simple  operator  such  that  div  curl  B  vanishes 
identically.  This  problem  is  solved  by  inclusion  of  curl  B 
in  Eq.  (4.4)  . 

4.2.2  Code  Structure 

A  block  diagram  of  the  hybrid  code  is  shown  in  Figure 

4.4.  Note  the  use  of  the  Incomplete  Cholesky  Conjugate 
r  q  i 

Gradient  (ICCG)  method  to  solve  large  sets  of  linear 
equations  in  both  major  code  modules.  It  is  the  availability 
of  this  highly  efficient  algorithm  that  makes  a  code  of  this 
type  possible. 

The  computational  grids  are  generated  using  the 
MESHGEN^~^  code.  (This  code  has  proved  to  be  less  than 
satisfactory.  We  expect  that  new  mesh  generators  soon  to 
be  available  will  be  substantial  improvements.)  The  initi¬ 
alization  phase  uses  the  finite  element  technique  to  construct 
the  Poisson's  equation  matrices,  and  converts  the  initial 
potential  distribution  to  an  initial  charge  distribution. 
Conducting  objects  and  boundary  nodes  may  be  floating  or  held 
at  fixed  potential. 

The  electromagnetic  and  electrostatic/particle  modules 
use  ICCG  and  otherwise  standard  techniques  to  solve  the  equa¬ 
tions  described  above.  Each  module  maintains  its  own  fast 


43 


Figure  4.3.  Centering  scheme  for  transverse  electromagnetic 
code.  Scalar  quantities  (potential  magnetic 
field  B^,  current  divergence  V* J)  are  centered 
at  solid  points,  while  vector'quantities  (trans¬ 
verse  electric  field  E,  current  density  J)  are 
centered  at  crosses. 
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Block  diagram 


random-access  file  in  which  mesh  information,  particle  in¬ 
formation,  field  information,  etc.  are  stored.  These  files 
may  be  saved  for  later  restarts  or  other  investigations. 

The  code  prints  out  electrostatic  potential  information 
periodically  as  requested.  Particle  information  is  stored 
on  a  disk  file  for  later  plotting  of  trajectories.  Surface 
current  information  is  both  printed  and  stored  for  later  pro¬ 
cessing. 

4.3  SIMULATION  OF  SGEMP  EXPERIMENTS 

We  used  this  newly-developed  code  to  simulate  experi¬ 
ments  of  the  type  performed  by  MRC.  We  chose  not  to  attempt 
to  reproduce  any  particular  shot  because  the  combined  effects 
of  poorly  defined  pulse  shape  and  emission  spectrum  and  the 
nonuni form  precharge  patterns  would  vitiate  any  agreement  or 
lack  thereof.  Rather,  we  simulated  generic  shots  to  study 
the  variation  of  response  with  electron  energy  and  current 
in  the  general  range  of  these  experiments,  and  to  develop 
simple  estimates  for  the  peak  response. 

Our  first  results  were  given  in  detail  in  our  quarterly 

report  of  August  1979.^*^  In  those  calculations  we  emitted 

2 

from  the  insulator  a  square  pulse  of  10  mA/cm  of  electrons 
-9 

for  20  x  10  sec,  with  an  exponential  spectrum  of  character¬ 
istic  energy  300  eV.  We  found  that  the  highest  surface  cur¬ 
rents  occurred  at  the  edge  of  the  emitting  insulator,  and 
that  the  peak  values  increased  from  ^1.8  A/m  without  precharge 
to  ■'''li  A/m  with  the  insulator  differentially  precharged  to 
-10  kV,  an  enhancement  of  ^6.  Similar  enhancements  were 
observed  elsewhere,  e.g.,  on  the  back  corner  the  response 
increased  from  "vO.3  A/m  to  ^1-5  A/m. 

While  we  were  pleased  with  the  agreement  between  our 
calculated  surface  currents  and  those  measured  by  MRC,  we 
were  concerned  that  we  predicted  blowoff  (i.e.,  final 


potential  of  test  body)  far  in  excess  of  that  measured.  We 
predicted  that  the  test  body  would  rise  to  nearly  10  kV, 
rather  than  the  1-2  kV  that  was  measured.  The  main  reason 
for  this  difference  was  that  our  estimate  of  the  capacitance 
of  the  test  body  to  the  tank  (^75  pf)  was  far  below  MRC's 
value  of  240  pf.  The  discrepancy  was  due  in  part  to  tank 
clutter  and  instrumentation,  in  part  to  our  neglect  of  rear 
surface  fields  in  the  electrostatic  grid  (Figure  4.1),  and 
in  part  because  neglect  of  potential  gradients  along  the 
surface  caused  excessive  blowoff. 

To  do  more  highly  resolved  calculations  we  chose  a 
baseline  case  as  follows: 

“8 

Pulse  shape:  square  wave  of  width  2  x  10  sec 

Voltage  profile:  V(r)  =  VQ  Jq  (2.4  r/R)  (r  <  R) 

2 

Emission  current  (J) :  200  A/m 

Electron  spectrum  (T) :  30  eV  exponential 

Here  J  is  the  ordinary  Bessel  function,  and  VQ  was  chosen  as 
-10  kV.  V  ,  J,  and  T  were  varied.  Zones  of  0.2  cm  immediately 
in  front  of  the  emitting  surface  were  found  to  resolve  the 
space  charge  layer  nicely.  Some  results  are  shown  in  Table 
4.1.  Appendix  B  contains  a  selection  of  surface  current 
response  curves  and  particle  trajectory  plots. 

The  main  conclusion  to  be  drawn  from  Table  4.1  is  that 
there  is  a  well-defined  limit  for  SGEMP  response  to  low 
energy  electrons  under  differentially  precharged  conditions. 

In  this  limit  the  response  depends  only  weakly  on  the  electron 
current  density  and  spectrum,  and  strongly  on  the  precharge 
potentials.  Note  also  that  both  the  surface  currents  and 
blowoff  are  in  good  agreement  with  MRC  results. 
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An  instructive  way  of  presenting  the  spatial  dependence 
of  the  skin  current  beneath  the  insulator  is  as  follows.  A 
rough  maximum  is  given  by  supposing  that  all  current  emitted 
within  radius  r  [nr^J]  is  replaced  by  surface  current  across 
r  [2irrK].  Thus  we  have 


K  (r) 
max 


(4.5 


The  ratio  K/Kmax  is  then  indicative  of  the  suppression  of 
response  due  to  space  charge  effects. 

Figure  4.5  shows  K(r)/'K  (r)  for  the  cases  run.  Note 

ITlclX 

the  modest  rise  of  the  baseline  curve  (I)  as  the  particle 
energy  is  increased  by  a  factor  of  16  (II,  III),  and  the 
sharp  drop  as  current  increases  by  only  a  factor  of  2  (IV) . 
Increasing  the  potential  (IV,  V)  markedly  increases  the  ef¬ 
ficiency  of  low  energy  electrons  in  producing  a  response. 

An  alternative  means  for  studying  the  surface  current 

spatial  dependence  is  in  terms  of  conductivity.  Under  this 

contract  we  developed  a  theory  for  the  conductivity  of  a 

space  charge  layer  (Appendix  A) .  In  this  case  there  is  a 

double  layer  of  space  charge,  divided  by  the  virtual  cathode 

which  forms  above  the  emitting  surface.  The  lower  layer  has 

-5  1/2 

a  conductivity  :  1  x  10  T  .  The  upper  is  characterized 

by  an  energy  V  J  (2.4  r/R) ,  and  a  lower  coefficient  due  to 

°  °  -6  1 
having  a  unidirectional  current:  c>2  -  5*10  [Vq  Jq(2.4  r/R)  ]  j 

The  surface  current  is  the  image  of  this  conductivity  current, 

and  is  thus  given  by 

K (r )  =  (c^  +  o2)  |VQ!  (2.4/R)  J1(2.4r/R)  .  (4.6) 


We  found  the  computed  peak  currents  could  be  approximated 
using  the  parameters 


/2 
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(4.7) 


oL  =  8.2  x  10  6  T1/2 

0o  =  3.27  x  IQ-6  [V  J  {2.4  r/R)]1/2  . 

2  o  o 

These  curves  are  shown  in  Figure  4.6.  They  reproduce  well  the 
peak  heights  and  their  variation  with  energy,  and  the  outward 
movement  of  the  peaks  with  increasing  energy.  They  predict  a 
peak  at  a  somewhat  smaller  radius  than  the  actual  computer 
simulations . 

4.4  SIMULATION  OF  SPACE  CHARGE-LIMITED  DISCHARGES 

The  SGEMP  code  described  in  the  previous  sections  can 
readily  be  modified  to  simulate  space  charge-limited  dis¬ 
charges.  To  do  this  ve  change  the  electron  emission  routine 
to  emit  that  current  cf  zero-kinetic-energy  electrons  which 
maintains  a  non-negative  surface-normal  electric  field.  This 
current  is  given  by  the  Child-Langmuir  law: 

J  =  (4co/9)  (2e/m} 1/2  IeJ  3/2/{Az)  1/2 

where  Lz  is  the  mesh  spacing  in  front  of  the  emitting  surface, 

and  E  is  the  mean  electric  field  across  &z  (assumed  electron 
z 

repelling} .  Not  only  is  this  treatment  applicable  to  discharges 
per  se,  but  it  is  equally  appropriate  to  cases  in  which  x-rays 
cause  a  plethora  of  low  energy  electrons  to  be  emitted  from  a 
precharged  surface. 

In  connection  with  a  theory  of  surface  discharges. 

Section  5.6  presents  some  emission  current  results  for  the 
MRC  geometry  assuming 


r  <m) 

4.6.  Predictions  of  peak  SGEMP  surface  currents  using 
Eqs .  (4. 6-4. 7)  (solid  curves).  Dot-dash  curves 
are  Cas's  IV  and  V  of  Table  4.1. 


where  V  =  -15  kV.  Here  we  concentrate  on  the  surface  cur- 

o 

rents  associated  with  such  a  discharge.  Parameters  character¬ 
izing  this  response  are  shown  in  Table  4.2.  For  the  Bessel 
function  case  the  response  is  similar  to  the  SGEMP  calcula¬ 
tion,  the  latter  giving  higher  surface  currents  due  to  non- 
negligible  effects  of  a  finite-energy  emission  spectrum.  As 
the  healing  length  parameter,  A,  becomes  short,  the  response 
rises  markedly.  Also,  a  new  effect  is  seen,  viz.,  a  negative 
precursor  in  the  surface  current  interior  to  the  dielectric 
(Figure  4.7).  This  can  be  explained  by  examining  Maxwell's 
equation 

VxH  =  J  +  e  |r-  E  . 

~  ~  ~  O  at  ~ 

In  this  interior  region  the  second  term  on  the  right,  caused 
by  increased  emission  at  larger  radii,  initially  dominates 
the  first,  causing  the  observed  effect. 
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SECTION  5 


A  THEORY  OF  DIELECTRIC  SURFACE  DISCHARGES* 

There  exists  a  large  body  of  experimental  data  which 
measure  peak  blowoff  currents  from  the  discharges  of  space¬ 
craft  dielectric  materials.  In  order  to  predict  how  large 
surfaces  behave,  several  researchers^12  14 ^  have  measured  how 
the  peak  blowoff  currents  scale  with  area.  Attempts  at  a 
self-consistent  theory  to  predict  the  observed  scaling  have 
all  suffered  from  qualitative  weaknesses ^1~>^  and  have  not 
been  successful  quantitatively.  In  this  section  we  present 
a  simple  model  of  surface  charge  emission  which  reproduces 
the  observed  area  dependencies  of  surface  discharge  quite 
well.  The  model  predicts  spatial  electron  currents  and 
therefore  is  useful  in  estimating  induced  surface  currents 
for  a  wide  variety  of  configurations.  We  make  no  attempt  to 
analyze  microscopic  phenomena  within  charged  dielectrics  or 
to  predict  discharge  thresholds.  While  limiting  the  overall 
scope  of  our  model,  the  independence  from  internal  dielectric 
microprocesses  allows  our  theory,  in  principle,  to  be  appli¬ 
cable  to  all  good  dielectrics,  given  only  surface  voltages 
prior  to  breakdown . 

5.1  REVIEW  OF  THE  DATA 

Several  researchers^12  1^'1^  have  exposed  thin  dielec¬ 
tric  samDles  to  electron  beams  which  caused  arcing.  We  shall 

[121 

focus  on  the  area  scaling  experiments  performed  by  Balmain. 

In  these  experiments  circular  Mylar  samples  were  exposed  to  a 
20  kV  electron  beam  in  a  high  vacuum  tank.  The  samples  were 
masked  by  grounded  aluminum  plates  to  minimize  edge  inhomo¬ 
geneities.  The  beam  density  was  of  the  order  of  1  pA  per 
square  centimeter.  Balmain  found  that  essentially  all  the 
* 

For  presentation  at  IEEE  1980  Nuclear  and  Space  Radiation  Ef- 
fects  Conference,  Ithaca,  NY,  July  1980. 
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surface  charge  was  blown  off  to  the  tank  wall  and  could  be 

measured  as  current  to  the  grounded  sample  backing  plate* 

The  peak  amplitude  of  this  current  scaled  approximately  as 

2 

the  area  for  small  samples,  A  <_  0.1  cm  ,  and  as  the  square 

2 

root  of  the  area  for  large  samples,  A  >  1  cm  .  The  tank 
walls  were  a  minimum  of  about  15  cm  from  the  sample.  Cur¬ 
rents  observed  ranged  from  0.01  amperes  to  100  amperes,  with 
pulse  durations  10-100  ns. 

5 . 2  THE  MODEL 

The  most  striking  feature  of  the  discharge  data  was 

the  magnitude  of  the  observed  currents.  For  the  1  square 

centimeter  sample  the  current  was  10  amperes.  Given  the 

distance  between  the  sample  and  the  tank  wall,  if  the  sample 

had  been  an  isolated  cathode.  Child's  law  (one-dimensional 

space  charge  limited  emission)  would  have  predicted  a  maxi- 
-2 

mum  of  3  x  10  amperes.  While  neutralization  of  the  electron 

space  charge  by  an  expanding  surface  plasma  could  allow  high 

6 

currents,  the  velocity  of  a  plasma  blowoff,  'vlO  cm/sec,  is 
far  too  slow  to  form  a  conducting  channel  within  the  observed 
pulse  duration  of  a  few  hundredths  to  a  tenth  of  a  micro¬ 
second.  In  order  to  achieve  th=  magnitude  of  the  currents 
observed,  it  is  apparent  that  most  of  the  sample  surface  must 
participate  in  the  discharge. 

The  model  we  propose  relies  on  the  edge  potential 
gradients  to  enhance  the  space  charge  limited  emission  from 
the  surface.  We  assume  that  some  observed  arc  discharge 
occurs  within  the  dielectric.  This  discharge  produces  enough 
radiation  to  trigger  space  charge  limited  emission  simul¬ 
taneously  from  the  entire  surface.  The  large  magnitude  of 
the  ejected  electron  current  results  from  the  strong  electric 
fields  near  the  edge  of  the  dielectric.  Below  we  present  a 
quantitative  formulation  which  predicts  both  small  and  large 
limits  and  the  transition  between  them.  The  model  is,  in 
principle,  entirely  predictive. 

,  ) 
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Before  dealing  with  an  actual  experimental  configura¬ 
tion,  it  is  instructive  to  examine  the  space  charge  limited 
current  emission  between  two  semi-infinite  conducting  half 
planes  (Figure  5.1).  This  problem  has  an  analytical  solution 
due  to  the  absence  of  a  scale  length  (see  Section  3) .  The 
first  result  is  that  the  density  of  emission  current  increases 
as  the  square  of  the  distance  to  the  edge 

V3/2 

J  -  (5.1) 

X 

This  corresponds  to  a  strongly  singular  current  near  the  edge. 
The  implication  for  experiments  is  that  if  the  voltage 
drops  off  at  the  edge  of  a  sample  on  a  scale  small  compared 
to  the  sample  radius,  then  the  emission  near  the  edge  will 
dominate  the  total  emission.  The  circumference  is  then  the 
important  emitting  region.  It  follows  that  the  current  is 
proportional  to  the  square  root  of  the  sample  area. 

A  second  important  result  for  the  semi-infinite  plane 
example  is  that  particles  travel  a  large  distance  parallel  to 
the  surface  before  being  collected.  A  particle  emitted  a 
distance,  d,  from  the  edge  of  the  cathode  lands  on  the  anode 
approximately  38  d  from  the  edge.  This  implies  that  a  small 
grounded  mask  will  collect  the  charge  emitted  only  from  a 
tiny  region  of  dielectric.  Thus,  most  of  the  charge  will 
travel  to  the  vacuum  tank  walls. 

The  obvious  drawback  of  the  step  function  potential 
calculation  is  the  singular  total  emission  current.  Certain 
other  semi-infinite  problems  without  scale  length  are  solvable 
which  yield  bounded  currents  per  unit  edge  length.  However, 
they  are  not  easily  related  to  realistic  cases.  In  Sections 
5. 3-5. 5  we  present  a  theoretical  formulation  for  emission  at 
an  edge  with  a  single  scale  length.  Numerical  methods  must 
be  used  to  predict  actual  current  profiles.  We  found  that 


straightforward  use  of  the  hybrid  SGEMP  code  described  in 
Section  4  was  the  best  tool  for  this  purpose.  Section  5.6 
contains  results  for  an  emitting  disk  on  a  ground  plane, 
which  is  a  good  model  for  comparison  with  experiment. 


5.3  GENERAL  DIMENSIONAL  CONSIDERATIONS  FOR  SPACE  CHARGE 
LIMITED  EMISSION 

We  consider  a  general  geometry  with  electrostatic 
potential 


V (r)  =  Vqv (r/X)  =  Vqv(z)  (5.2) 


where  V(r)  is  specified  on  the  boundaries,  some  subset  of 
which  are  capable  of  emitting  electrons  in  space  charge 
limited  fashion,  Vq  characterizes  the  magnitude  of  the  poten¬ 
tial,  and  X  characterizes  the  spatial  variation  of  potential 
on  the  boundaries.  The  space  charge  density  is  given  by 

P(r)  =  -eo?2V(r)  =  -  (e^/X2)  V2v (z)  (5.3) 

and  the  current  density  by 


J(r) 


p(r)U(r)  =  -U0V0/X2) 


y3/2 


o 


j  (z) 


(5.4) 


Here  we  have  assumed  non-crossing  particle  orbits  so  that 
U(r)  is  a  unique  function  of  space,  and  we  introduce 


(5.5) 
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and 


z  = 


dz 

dr 


(5.6) 


is  the  dimensionless  particle  velocity. 

The  above  allows  the  entire  problem  to  be  reduced  to 
dimensionless  equations: 

The  steady-state  requirement 


V  •  j  =  0 


(5.7a) 


Poisson's  equation 


2  # 

V^V  =  +jjj/jz|  (assuming  electron  currents)  (5.7b) 


and  the  equation  of  motion 


z  =  2  • 


(5.7c) 


These  three  equations,  together  with  the  boundary  conditions 
and  the  scaling  relations 


V  (r) 


=  Vqv (r/X ) 


(5.8a) 


2  v1/2  ..3/2 


/  2es  \ 

iri 


j (r/X) 


(5.8b) 


e  v 

p(r)  =  -  -2-0  ! j (r/X) |/|z(r/>.) 
X^  -  “ 


(5.8c) 


then  permit  a  complete  (numerical)  solution  of  the  problem. 
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5.4  PROBLEMS  WITH  NO  LENGTH  SCALE 

We  consider  here  a  general  class  of  problems  relating  to 
the  breakdown  near  an  edge.  Suppose  two  semi-infinite  planes 
intersect  at  angle  0q  (Figure  5.2)  and  have  potentials 

V(?)  =  Vo(|ri/X)V  fv(0)  {5  o a. 


fv(0)  =  0 


(5.9b) 


This  problem  can  be  seen  to  have  no  length  scale  because  it 
is  invariant  to  the  transformation 


A  -*■  <x  A 


V  -  a  v 
o  o 


Using  this  in  (5.8b)  gives 


J(r) 


■(5s) 


1/2 


|v  ! 

1  o' 


3/2 

3v/2-2  . 

—  a  /  j  (r/aA ) 


or 


J(r,6)  = 


■m 


2,1/2  |V  |V2 
1  o' 


x  2 


l*f/2-2  jci.ej 


(5.10a) 

(5.10b) 


(5.11a) 


(5.11b) 


or 


J(rf3)  = 


•(sf(ii) 


3/2 


,3v/2-2 


3(1,9) 


( 5 . 11c) 
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Table  5.1.  Emission  Current  Density  (Relative  to  Child- 

Langmuir-like  Value)  and  Anode  Current  Density 
Reduction  Factor  for  Geometry  of  Figure  5.2  with 
Step-function  Potential. 


0 

o 

9<V 

CD 

o 

0.1 

1.0 

1.003 

0.5 

1.0 

1.08 

1.0 

0.95 

1.37 

7T/2 

0.90 

2.29 

2.0 

0.81 

4.16 

2.5 

0.63 

10.5 

3.0 

0.45 

28.9 

TT 

0.38 

38.2 

0Q  =  angle  between  half  planes. 

g(0Q)  =  dimensionless  parameter  from  Eq.  (5.12) 
used  to  calculate  -pace  charge  limited 
current  between  infinite  half  planes. 

R(0q)  =  distance  from  origin  of  electrons  striking 
anode  when  emitted  from  cathode  at  R  =  1.0. 


5.5 


A  CLASS  OF  PROBLEMS  WITH  ONE  LENGTH  SCALE 


We  consider  a  class  of  problems  similar  to  the  above, 
but  require 


0  =  TT 

O 


V  <  0 
O 


V(r,o)  =  0 


V(r,ir)  =  Vq(1  -  e 


or 


=  Vq  tanh  (r/A) 


(5 

(5 

(5 

(5 


(5 


or 


2V 
_ o 

IT 


arctan 


(5 


or  some  similar  function.  For  a  point  on  the  cathode 
r  >>  A,  any  problem  of  this  class  looks  like  the  v  =  o 

(step  function)  problem  of  the  preceding  section,  and  thus 

_2 

has  emission  current  J  «  r  .  On  the  other  hand,  for  r  <<  A 

they  all  reduce  to  the  v  =  i  problem,  and  thus  emit  zero 

current.  J(r)  thus  has  a  maximum  for  r  -v  A .  The  value  at 

this  maximum  can  be  far  in  excess  of  the  value  g(n)  v2/^2  A-2 

o 

inferred  from  (5.12),  since  the  current  is  limited  not  by 
its  own  space  charge  (which  is  rapidly  removed  by  the  tan¬ 
gential  electric  field)  but  by  space  charge  from  the  more 
weakly  emitting  region  r  >  A . 


.14a) 

.14b) 

.14c) 
.  14d) 

.  14e) 

.  14f ) 


5.6 


A  PROBLEM  WITH  TWO  LENGTH  SCALES 


In  realistic  problems  there  tend  to  be  two  scales  of 
length:  one  characterizing  the  geometry  and  one  characteriz¬ 
ing  the  surface  voltage  profile  near  ar  edge.  We  consider 
in  particular  a  charged  dielectric  disk  of  radius  R  on  a 
ground  plane.  It  is  apparent  the  Eqs.  (5.8a,b)  must  now 
be  written 

V (r)  =  VQv(r/X;  R/X) 

/2e£2\V2  V3/2 

J<E>  =\~ T2/  j(r/X;  R/X) 


(5.15a) 


(5.15b) 


We  can  simplify  (5.19b)  by  integrating  over  the  disk  surface: 


R 

Itot  -/  2”rdrJ<r» 
0 


R/X 

j  (z;R/X)  zdz 

0 

(5.16) 


Written  thus,  it  is  manifestly  clear  that  the  total  current 
3 /2 

is  VQ  times  a  function  of  R/X  only. 

We  have  used  the  self-consistent  electromagnetic  parti¬ 
cle  code  described  in  Section  A  to  simulate  this  problem  for 
a  disk  with  R  =  40  cm  and  a  potential  of  the  form 


V  (r) 


V 

o 


1  -  e~(R”r)/X 
1  -  e"R/X 


(5.17) 


This  potential  has  a  peak  value  VQ  at  the  center  of  a  disk  of 
radius  R,  and  goes  to  zero  at  the  edge  with  a  scale  length, 

X. 
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dV  .  o  1 

3F  <r=R)  "  —  - - =57x7  (5-13> 

(1  -  e  '  ) 

/  VQ/X  R  >>  X 

(  Vq/R  R  <<  X  . 


The  parameter  X  was  varied  from  3  cm  to  infinity  (triangle 
potential).  The  results  are  shown  in  Figures  5. 3-5. 6.  For 
X  <<  R,  it  is  apparent  from  Figure  5.3  that  the  maximum  cur¬ 
rent  density  is  given  by 

J  (X<<R)  Z  2.45  x  10“7  V3/2/X2  (5.19) 

max  o 

in  keeping  with  Eq.  (5.8b),  and  from  Figure  5.6  that 

Itot  s  +  2.9X10'6  f)  V3/2  (5.20) 

in  keeping  with  Eq.  (5.16).  The  current  maximum  occurs  at 
rmax  ~  R  ”  1.3X,  and  the  emission  vanishes  at  r  '  R  -  X/2 
in  keeping  with  the  ideas  of  subsection  5.5. 

5.7  COMPARISON  WITH  EXPERIMENT 

A  simple  comparison  with  experiment  is  presently  not 
possible  since  the  sample  surface  voltage  profile  just  prior 
to  discharge  was  not  measured.  We  would  hope  that  voltage 
profiles  would  become  available  soon.  However,  if  we  assume 
the  analytical  form  of  the  potential  in  Eq.  (5.17),  then  for 
any  example  there  are  only  two  physical  parameters  to  be 
estimated:  VQ  and  X.  Balmain's  data  extrapolate  to  10  amps 
for  aim2  area  sample.  If  we  assume  VQ  =  -15  kV,  Eq.  (5.20) 
gives  X  =  0.3  cm.  To  extrapolate  to  small  samples  we  suppose 
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Figure  5.6.  Total  space  charge  limited  emission  current  from 
disk  having  potential  profile  Eq.  (5.2)  with 
V  =  -15  kV. 


dV 

dr 


r=R 


constant 


X  =  constant 


(5.21) 


at  breakdovm.  Using  the  above  data  we  find  dV/dr|r_R  = 

5  x  10®  V/m.  The  predictions  of  Eq.  (5.17),  (5.20),  and 

(5.21)  for  the  current,  together  with  Balmain's  peak  current 
T 121 

data  are  shown  in  Figure  5.7.  We  suggest  that  the  dis¬ 
crepancies  for  the  very  small  examples  are  due  in  part  to  the 
short  pulse  length,  which  was  comparable  to  the  inverse  band¬ 
width  of  Balmain's  experiment.  For  large  samples  R/X  >>  1, 

1/2 

Eq.  (5.22)  predicts  I  «  R/X  ■«  A  '  .  For  small  samples, 

Eq.  (5.3)  gives  Vq  =  R  dV/dr,  so  that 


(If  the  blowcff  is  considered  emission  limited,  we  would 
have  I  “  A. ) 

While  this  comparison  involves  fitting  a  parameter  and 
assumed  potential  form,  the  asymptotic  behavior  and  the  magni¬ 
tude  agreement  lend  credence  to  the  theory.  Definitive  tests 
performed  with  measured  and  well  characterized  voltage  pro¬ 
files  would  be  a  tremendous  aid  in  evaluating  the  hypotheses. 


5.8  CONCLUSION 

We  have  presented  a  model  fo i.  the  "blowoff"  of  surface 
charge  from  a  dielectric  surface  during  discharge.  This  model 
successfully  explains  both  the  magnitude  of  charge  release 
currents  and  how  such  currents  scale  with  sample  area.  While 
the  specific  dielectric  mechanisms  which  trigger  discharges 
are  not  included,  their  absence  makes  the  results  generaliz- 
able  to  a  wide  range  of  materials - 
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Emission  current  calculated  using  assumption  (5.21)  and  Eq.  (5.20),  com¬ 
pared  with  data  of  Balmain  (Reference  2)  . 


One  feature  of  the  model  is  that  the  entire  area  simul¬ 
taneously  discharges.  This  accounts  for  the  substantial  loss 
of  surface  charge.  The  model  accurately  predicts  the  circum¬ 
ference  dependence  of  discharge  currents  for  moderately  large 
samples.  Extensions  of  the  model  to  include  finite  return 
path  impedance  should  account  for  observed  impedance  effects. 

Finally,  since  emitted  electron  trajectories  are  cal¬ 
culated,  the  model  predicts  the  effective  electromagnetic 
driver  of  surface  discharges.  Using  such  a  driver  and  SGEMP 
codes,  one  can  (Section  4.£)  predict  the  response  of  space¬ 
craft  to  electrostatic  surface  discharges. 
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A  mechanism  for  simultaneous  flashover 
of  an  array  of  dielectric  surfaces  is  pro¬ 
posed.  The  timescale  and  current  level  of  the 
resultant  discharge  are  evaluated,  and  conse¬ 
quences  for  SGEMP  response  are  discussed. 

Introduction 

In  this  paper  we  examine  how  the  dis¬ 
charge  of  thin  dielectric  surfaces  can  sig¬ 
nificantly  enhance  SGEMP  response.  The 
mechanism  of  discharge  involved,  viz.  surface 
flashover,  has  not  previously  been  discussed 
in  the  SGEMP  context,  although  apparent 
SGEMP-triggered  discharges  were  observed  in 
tests  on  the  Skynet  qualification  model.1 

The  basic  problems  of  SGEMP  are 
generation  of  electromagnetic  energy  in  re¬ 
sponse  to  an  x-ray  pulse  and  coupling  of  this 
energy  into  the  structure  and  circuitry  of  a 
spacecraft.  Most  studies  have  focused  on 
electromagnetic  signals  produced  directly  by 
2  4 

the  10  -10  eV  primary  photoemitted  elec¬ 
trons.  Such  signals  are  not  independent  of 
the  initial  state  of  the  test  object.  Uni¬ 
form  negative  precharging  leads  to  a  slightly 
enhanced  response.2  Substantial  further  en¬ 
hancement  due  to  differential  precharging 
have  previously  been  examined3  and  effects  of 
this  type  are  currently  under  study. 

Yet  larger  responses  would  result  from 
release  of  stored  electrical  energy  as  a 
synergistic  effect.  Many  common  insulating 
spacecraft  surfaces  have  capacitances  c  to 
2 

ground  of  n-10  pf/cm  .  Differential  charging 
of  such  a  surface  to  o  =  10  kv  thus  produces 

12  -3  2 

an  energy  density  t  =  ^  C;  =  5  x  10  J/cm  . 

For  an  incident  x-ray  fluence  of  this  amount 
the  electromagnetic  energy  den  :ity  that  can 
be  created  by  primary  photoelectrons  is  about 
a  thousand  times  smaller  (even  less  if  the 
emission  current  is  limited  by  space  charge) . 
Thus  a  discharge  of  insulating  surface  layers 
can  lead  to  a  substantial  response  even  if  it 
couples  very  inefficiently  to  electromagnetic 
energy. 

The  most  rapid  discharge  of  a  substantial 
dielectric  area  would  occur  by  surface  flash- 
over.  By  this,  we  mean  that  an  electron  layer 
forms  just  above  the  surface,  causing  it  to 
become  effectively  conducting,  i.e  ,  capable 
of  carrying  large  currents  parallel  to  the 
surface.  In  this  paper  we  address  the  ques¬ 
tions:  (1)  If  such  a  layer  is  formed,  what 
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are  the  characteristics  of  the  subsequent 
discharge?  (2)  Under  what  conditions  might 
such  a  discharge  be  initiated  by  an  x-ray 
pulse? 

Surface  Conductivity 

Let  us  begin  by  addressing  the  question 
of  the  surface  conductivity  associated  with 
a  space  charge  layer.  If  an  emitted  electron 
returns  to  the  surface  after  a  time  t ,  it 
will  travel  along  the  surface  with  mean  ve- 
1  e 

locity  v  =  —  Em  t  ,  where  E„  (assumed  uni¬ 

form)  is  the  electric  field  parallel  to  the 
surface.  Such  a  field  results  from  non-uni¬ 
formity  of  the  precharging  due,  for  example, 
to  ambient  surface  conductivity  near  edges4 
(see  Figure  1) .  The  area  density  of  electrons 
above  the  surface  emitted  with  energy  between 

c  and  c  +  dc  is  t-j(c)  i  (e)de .  f  .is  a  surface 
d  C 

current 


(i)  - 


||(E)T2(e)  dc 


is  obtained,  or,  a  surface  conductivity 


c  (ohms 


h--*-  f  i 
“  J.  1 


!^-(e)t2-(c)  de.  (2) 


DIELECTRIC— i 


I—  DIELECTRIC 


Figure  1.  Schematic  diagram  of  potential 
contours  near  edges  of  charged  dielectrics. 
Reasonable  value  for  potential  gradient  along 

dielectric  surface  might  be  1  kv/l  mm,  or  10 

V/m. 


we  can  evaluate  this  conductivity  for 
either  the  monoenergetic  (Child-Langmuir) 
case  j'U)  "v  6  (c— c^),  or  for  the  thermal  case 

j'(c)  ‘v  e_C//T.  (We  neglect  all  initial 
transverse  velocities.) 
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For  the  monoenergetic  case  (see  Appendix 
A)  the  flight  time  t  is  given  by 


t  (sec)  =  1.09  x  10-8  el/4  J-1/2 
and  from  Eq.  (2)  we  obtain 

a (ohms'1)  =  1.05  x  10~S  e1/2 


(3a) 


(3b) 


where  e  is  emission  energy  (eV)  and  J  is 
2 

emission  current  (A/m  ) .  Note  that  o  is  in¬ 
dependent  of  J  (although  the  conductivity 
will  saturate  if  J  is  too  small).  Thus  the 
surface  resistivity  of  a  monoenergetic  space 

charge  layer  is  ^95  kO/e1^2. 

For  the  thermal  case  we  find  (see 
Appendix  B) 


1/2 

x(e)  =  (g*)1/4  C-f)  exp (e/2T 


x  erf[(e/2T)1/2] 


=  7.27  x  io"9  t1/4  J_1/2  exp(e/2T) 


(4a) 


■/ 


p  dz 


(7) 


and  to  integrate  the  current  density  j  over 
the  vertical  extent  of  the  space  charge  layer 


K„  =y  j  dz  . 


This  allows  us  to  write 


U  +  Y 


K„  =  °  - 


(8) 


(9) 


The  areal  charge  density  can  be  considered  as 
having  two  components 


q  =  qs  +  qv 


(10) 


'  here  in  the  charge  residing  at  the  surface 
and  q  is  the  contribution  in  the  space  charge 

limited  layer.  Under  conditions  of  SGEMP 
relevance  the  space  charge  limited  current 

density  is  typically  J  <  10  2  A/cm2  and  cor¬ 
responding  flight  times  7  10~10  secs  so  that 


x  erf  ( (e/2T) 1/2] 


(4b) 


q  \  J  T 


where  T  is  the  temperature  (eV)  and  CQ  the 
permittivity  of  free  space. 

To  evaluate  tha  conductivity  we  must  in¬ 
sert  into  Eq.  (2)  a  convergence  factor 
exp(-e/Tc),  reflecting  the  fact  that  very 

high  energy  electrons  leave  the  influence  of 
the  transverse  field.  We  then  find 


Ofohms'1)  =  4.6  x  10'6  T1/2  ( 


-12  2 
10  coulombs/cm 


By  comparison,  the  surface  charge,  qs,  is  much 

larger  since,  for  a  dielectric  layer  of  thick- 
-2 

ness  d  t  10  cm  the  capacitance  is 


C  = 


9  x  io 


-14 


■v  10  **  farads/cm2 


Thus,  the  surface  charge  density  is 


4  /  T  \1/2  /  T  \1/2 

n  \2T+T  /  arctan  \2t+t  / 


(5) 


(The  quantity  in  square  brackets  becomes 

unity  for  T  >>  T. ) 
c 

Surf  ice  Diffusion  Discharge 

We  turn  now  to  the  relationship  between 
conductivity  and  discharge  characteristics. 
The  space  charges  and  currents  must  satisfy 
the  equation 


! 


j  =  0 


(6) 


Since  there  is  negligible  net  current  in  the 
z-direction,  the  current  divergence  is  en¬ 
tirely  lateral  -  that  is,  parallel  to  the 
emitting  surface.  It  is  thus  appropriate  to 
define  an  areal  charge  density 


%  =  c* 


-7  2 

■v  10  coulombs/cm 


(111 


if  the  surface  is  charged  to  10  kV  with  re¬ 
spect  to  the  underlying  conducting  surface. 
Inserting  this  expression  for  q  and  the  rela¬ 
tion  implied  by  Eqs.  (1)  and  (2),  namely. 


into  Eq.  (9) ,  we  obtain 


C  ||  =  V  •  (oV$) 


£  |4  =  72*  . 

o  3t 


(12) 


(13) 
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(We  neglect  spatial  dependence  of  the  con¬ 
ductivity,  which,  in  the  treatment  used  nere, 
is  a  function  only  of  the  electron  emission 
spectrum.)  The  quantity  C  has  the  dimension 

2 

of  capacitance  per  unit  area  (farads/cm  )  and 

a  those  of  conductance  (ohms  so  that  the 
inverse  diffusion  coefficient  C/o  has  dimen¬ 
sions  of  time  per  unit  area.  The  discharge 
time  is  then  proportional  to  the  sample  area 
and  inversely  proportional  to  its  thickness. 

The  solution  to  (13)  for  a  1  cm  radius 
-7  2 

disk  with  C/o  =  5  «  10  sec/cm  is  shown  in 
Figure  2.  These  parameters  are  appropriate 
to  a  charged  solar  cell  cover  with  a  ground 
electrode  along  its  edge.  It  is  seen  that 
the  surface  is  about  20  percent  discharged  in 
—8 

10  seconds.  Simultaneous  discharge  of  many 
such  solar  covers  could  produce  a  strong 
electromagnetic  response. 


u  3x10 

m 

<A 


-V=9  kV 
-V=8  kV 
-V=7  kV 


.8  1.0 


R  (cm) 


Figure  2.  Contours  of  constant  potential  as 
functions  of  position  and  time. 


Secondary  Electron  Avalanche 

A  discrepancy  in  the  above  mechanism  ap¬ 
pears  when  we  examine  the  current  level  re- 

2 

quired  for  such  a  discharge.  If  a  10  pf/cm 
surface  at  10  kV  is  to  discharge  20  percent 

_  o 

in  l*  sec,  a  net  emission  current  density 
2 

of  2  A/cm  is  required.  Since  SGEMP  photo¬ 
electron  currents  are  much  less  than  this, 
some  current  amplification  mechanism  is 
necessary.  For  such  a  mechanism  we  suggest 
a  secondary  electron  avalanche,  which  was 
studied  in  a  related  context  by  Anderson.5 
The  basic  idea  is  that  electrons  striking  a 
surface  with  energy  in  the  range  ^50-1000  eV 
can  liberate  as  many  as  ■'-20  low  energy  (n*? 
eV)  secondary  electrons.  These  secondaries 
can  then  be  accelerated  by  the  electric  field 
parallel  to  the  surface  to  continue  the 
cascade. 

To  examine  the  plausibility  of  this 
mechanism,  we  return  to  Eq.  (4b).  Suppose  an 

2 

x-ray  pulse  produces  electrons  with  J  =  10 
2 

A/m  and  T  =  300  eV.  (We  do  not  mean  to  im¬ 
ply  here  that  a  thermal  cistribution 


adequately  represents  an  SGEMP  electron  emis¬ 
sion  spectrum.  Our  intent  is  simply  to  show 
that  the  high  energy  components  of  the  spec¬ 
trum  space-charge-limit  the  low  energy  elec¬ 
trons  in  such  a  way  thai  they  may  cause  a 
secondary-electron  avalanche.)  Since  erf (x) 

;  x  for  x  <  1,  the  flight  time  of  a  2  eV 
electron  is 


t  -v  7  x  10-9  (2  eV)  1/2  (300  eV)-1/4 


2  -1/2 
•  (100  A/m)  (2) 


-v  1.7  x  10  sec. 


In  a  lateral  field  of  3  x  io3  V/m,  it  will 

2  —4 

hop  a  distance  (e/2m)E„T  =  8  *  10  m,  and 

gain  240  eV  of  energy.  (This  field  is  ap¬ 
propriate  near  the  edge  of  a  small  sample 
charged  to  1  kV.)  It  will  then,  upon  landing, 
produce  "v-2-20  secondary  electrons.  In  fact, 
the  x-ray  produced  current  of  low  energy 
electrons  is  a  substantial  fraction  of  the 
total  current,  so  that  the  emission  current 
increases  considerably  with  just  one  hop. 

The  amplification  will  be  reduced  as  the 
secondary  electrons  become  limited  by  their 
own  space  charge.  The  steady-state  emission 
current  density  is 


,-5  1/2 


where  e  is  the  emission  energy  and  es  is  the 
energy  for  unit  secondary  production.  For 

e  =  2  eV,  e  =  30  eV,  E„  =  3  x  io3  V/cm  we 
s  2 

find  J  "v  4  A/cm  ,  far  in  excess  of  what  is 
ss  ' 

needed  to  sustain  the  discharge. 

Conclusions 

In  conclusion,  x-rays  can  trigger  the 
simultaneous  discharge  of  an  array  of  small 
dielectric  surfaces  through  the  mechanism  of 
secondary  electron  avalanche.  The  lateral 
electric  field,  resulting  from  the  difference 
in  potential  between  charged  dielectric  sur¬ 
faces  and  an  underlying  conductor  (e.g.,  as 
in  solar  panels  -  see  Figure  1),  provides  the 
energy  for  the  discharge.  The  SGEMP  effect 
(driven  by  the  displacement  current)  of  these 
discharges  will  be  greatest  when  the  segments 
2 

are  at  most  a  few  cm  in  area,  and  are  dif¬ 
ferentially  charged  to  1-10  kV.  The  dis¬ 
placement  current,  3E/3t  will  be  proportional 
to  the  dielectric  thickness.  By  way  of  com¬ 
parison  to  SGEMP,  the  currents  generated  by 
such  a  breakdown  are  (see  Eq.  (3b)) 


K„  =  oE„  -v  1  x  io"5  x  io6 


10  A/m 


for  conditions  where  an  electric  field  E  t 

10®  V/m  exists.  This  is  an  order  of  magni¬ 
tude  larger  than  the  SGEMP  surface  current 
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response  observed  for  fluence  conditions 

-1  -5  2  fi 

ranging  from  10  to  10  joules/cm  . 
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Appendix  A 

Conductivity  of  Monoenergetic  Space  Charge 
Layer 

The  potential  due  to  the  space  charge  of 
electrons  emitted  normally  and  monoenergeti- 
cally  (energy  c  (eV) ]  satisfies  the  equation 


J 

o 


-Be 


(B.l) 


where  e  =  1/2  mv2  is  the  normally  directed 

energy  component.  We  wish  to  calculate  (1) 
the  potential  4(z),  subject  to  boundary  con¬ 
ditions  4(0)  =  4'(“)  =  0;  the  flight  time 
x(e);  and  (3)  the  transverse  conductivity. 

The  space  charge  density  at  z  is  com¬ 
posed  of  contributions  from  those  electrons 
with  sufficient  energy  to  pass  z: 


CO 

p(z)  =  -2  J  ( (2/m)  (£+e4(z))l-1/2  de. 


-e6 (z) 


(B.  2) 


Using  (B.l)  and  integrating, 

o(z)  =  -  Jo(2nmB)1/2  eSe*  .  (B.3) 

The  potential  is  then  the  solution  of 


d2»i>  _  2J 
dz2  CQ/2e4/m 


(A.  1) 


o/eo 


J 

o 

e 

o 


(2nmB) 


1/2 


Bep 


(B.  4) 


A  virtual  cathode  forms  at  a  distance  d  from 
the  surface,  and  (A.l)  is  subject  to  boundary 
conditions  4(0)  -»  +e_,  4  (d)  =  0,  4'(d)  =  0. 
Letting  k  =  2/e  /2e7m,  we  find 


4 


(A. 2) 


subject  to  the  above  boundary  conditions. 
This  solution  is 


where 


d  =  1  e3/4  (kJ)-1/2 


The  flight  time  is  then  given  by 


(B.6) 


a 

2  f 

dz 

/  2  \1/4 

=  ei/«  j-1/2  (32m  co) 

J 

0 

Vir  *(z) 

\  e  / 

(A. 4) 

1.09  * 

10-8  c1/4 

J-1/2  . 

transverse  conductivity  is  then 

Inversely, 

z(4)/zq  =  /2  [e'6e$/2  -  1]  . 

The  flight  time  t  is  given  by 
z : (-e/e) 


(B.  7) 


(e)  =  2  J  1  (2/m)  (e+ep  (z) )  J-1/2  dz(B.8) 
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t  (e) 


e6e/2  erf  (0) 


Eeff(e)  “  E«  e 


(B.  12 


wr 


(B.  9) 


t  (e) 


=  7 


,27  x  lo"9  T1/4  J~1/2  e"e'/2Terf 
o 


c 


(A  better  defined  approximation  would  be  to 
treat  explicitly  the  dependence  of  trans¬ 
verse  field  on  height.  We  reserve  such  a 
treatment  for  future  work.)  Equations  (B.l) 
(B.9),  (B.ll),  and  (B.12)  then  give 


(B. 10) 

where  t  is  in  seconds,  T  in  eV,  and  J  in 

2  0 
A/m  . 


The  transverse  current  is  given  by 


o  =  e  (it/2m) 
o 


1/2 


CD 

/ 


-e/T 


c  erf2 


or 


(B.ll 


Kn  =  | J  /  |f  t2U)  Sef£(c)  dc  .  (B.ll) 
0 

Here  we  must  take  into  account  the  fact  that 
very  high  energy  electrons  leave  the  region 
of  the  transverse  field,  and  should  be  given 
reduced  weight  in  calculating  transverse  cur¬ 
rent.  To  achieve  a  simple,  one-parameter 
model  we  take 


(B.  14 
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RESULTS  FOR  SGEMP 

(BESSEL  FUNCTION  PROFILE)  CALCULATIONS 

We  present  here  particle  trajectory  plots  and  surface 
current  plots  for  the  five  SGEMP  cases  of  Table  4.1.  The 
particle  trajectory  plots  are  "open  shutter"  pictures  of 
electron  tracks.  Note  the  substantial  amount  of  charge  ac¬ 
celerated  radially  outward  by  the  surface  precharge  fields. 
There  is  also  a  substantial  current  layer  of  space  charge- 
limited  electrons  along  the  front  surface  extending  beyond 
the  edge  of  the  dielectric.  No  electrons  were  seen  to  hit 
the  side  or  rear  of  the  test  body,  though  a  few  may  do  so 
later  in  time. 

Surface  current  plots  are  presented  from  four  repre¬ 
sentative  surface  nodes  (see  Figure  B.l) .  Node  398  has  the 
peak  response  for  this  voltage  profile.  Node  400  is  at  the 
edge  of  the  dielectric.  Node  403  at  the  front  corner  and 
node  314  at  the  rear  corner  give  a  measure  of  ,he  electro- 
magnetically  propagating  response  likely  to  be  picked  up  by 
antennas  or  to  penetrate  through  apertures.  (The  surface 
currents  have  been  filtered  to  remove  high  frequency  numerical 
noise. ) 
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Figure  B.l.  Location  of  B-field  (surface  current)  node  points 
in  computer  model  of  differentially  precharged 
SGEMP  experiment. 
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Surface  currents,  Case  III 


Surface  currents.  Case  III 


Figure  B.8.  Particle  trajectories,  Case  IV 
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TECHNICAL  NOTE 

EFFECT  OF  EARTH'S  MAGNETIC  FIELD  ON 
PRECHARGING  EXPERIMENTS 

G.  W.  Schnuelle  and  M.  J.  Mandell 

The  ambient  magnetic  field  can  profoundly  affect  the 
charging  of  an  object  by  an  electron  gun  in  a  30  m  tank.  The 
Larmor  radius  for  electrons  is  rL  1.2  /e/B  meters,  where 
E  is  the  electron  energy  in  kV  and  B  the  magnetic  field  in 
gauss:  rL  ^  7  meters  for  a  10  kV  gun  with  B  =  0.5  gauss. 
Obviously  some  reduction  of  the  ambient  field  is  necessary. 

We  have  used  the  NASCAP  code  to  illustrate  the  effects  of 
various  levels  of  ambient  field  bucking. 

Figure  1  illustrates  the  object  used  in  the  simula¬ 
tions;  the  zone  size  is  1  meter.  Figures  2  through  4  show 
the  trajectories  of  a  1  milliamp  beam  of  10  kV  electrons  as 
they  travel  from  the  source  toward  the  object.  The  source 
is  positioned  16  meters  from  the  object,  and  the  magnetic 
field  is  oriented  out  of  the  plane  of  the  paper.  The  cal¬ 
culation  includes  the  combined  effects  of  the  ambient  mag¬ 
netic  field  and  the  electrostatic  field  from  the  charged 
object  on  the  trajectories.  In  Figures  2  through  4,  the 
object  is  uniformly  charged  to  -5  kV,  and  the  magnetic  field 
is  0.,  .1,  and  .25  gauss  (about  half  ambient)  respectively. 

With  the  field  bucked  to  half  ambient  (Figure  4) , 
electrons  cannot  reach  the  top  half  of  the  object.  With 
the  field  bucked  to  one-fifth  ambient  (Figure  3) ,  the  gun 
can  spray  the  entire  object,  h.  ,/ever,  even  in  this  ease 
the  field  causes  trajectories  reaching  the  top  of  the  object 
to  have  much  higher  curvature  than  those  striking  the  bottom. 
(Note  that  the  curvature  increases  as  the  electron  energy 
decreases  approaching  the  charged  object.)  As  a  result  of 
the  angle  dependence  of  secondary  emission,  the  net  current 
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density  on  the  top  portion  is  +7  x  10  A/m  ,  while  the 

“82 

current  at  the  bottom  is  -.4  x  10  A/m  .  Several  thousand 
volts  of  differential  charging  will  result  from  the  effects 
of  the  non-zero  magnetic  field. 

These  calculations  indicate  that  reduction  of  the 
ambient  field  to  approximately  0.1  gauss  should  suffice  to 
allow  10  keV  electrons  to  reach  a  test  object.  However, 
even  in  this  case  further  simulations  will  be  required  to 
estimate  the  effects  of  the  remaining  magnetic  fields  on 
differential  charging. 

To  reduce  the  horizontal  and  vertical  magnetic  field 
components  by  a  factor  of  10  over  a  large  tank  is  a  major 
task.  It  has  not  been  attempted  at  the  NASA  Johnson  Space 
Center;  instead,  experiments  have  been  performed  at  several 
field  strengths  and  the  results  extrapolated. 

Care  in  the  positioning  of  electron  guns  combined 
with  the  capability  of  moving  guns  on  pods  out  from  walls 
could  overcome  at  much  lower  cost  the  problem  of  precharging 
in  the  earth's  magnetic  field. 

An  alternative  method  of  achieving  differential 
charging  would  be:  (1)  bias  spacecraft  ground  to  +10  kV; 

(2)  turn  on  a  low  energy  plasma  source,  allowing  charge  to 
accumulate  on  dielectric  surfaces  until  they  are  near  tank 
ground;  (3)  return  spacecraft  to  tank  ground,  so  dielectric 
surfaces  are  at  -10  kV;  (4)  if  desired,  discharge  selected 
areas  with  UV  source. 
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PARTICLE  TRAJECTORIES 


Figure  3.  Trajectories  of  10  kV  electrons  traveling  16  m  to 
test  object.  Magnetic  field  is  0.1  gauss  (out  of 
plane  of  paper) .  Object  potential  =  -5  kV. 
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ATTN:  ATC-0 

BMD  Systems  Command 
Department  of  the  Army 
ATTN:  8DMSC-H 

Deputy  Chief  of  Staff  for  Rsch  Dev  &  Acq 
Department  of  the  Army 
ATTN:  DAMA-CSS-N 

Electronics  Tech  &  Devices  Lab 
U.S.  Army  Electronics  R5D  Command 
ATTN:  DRSEL 


DEPARTMENT  OF  THE  ARMY  (Continued) 

Harry  Diamond  Labs 
Department  of  the  Army 

ATTN:  DELHD-N-RBC,  R.  Gilbert 
ATTN:  DELHD-I-TL 

U.S.  Army  Cormiunications  Sys  Agency 
ATTN:  CCM- AD-LB 

U.S.  Army  Foreign  Science  &  Tech  Ctr 
ATTN:  DRXST-IS-1 

U.S.  Army  Missile  R&D  Command 
ATTN:  RSIC 

DEPARTMENT  OF  THE  NAVY 

Naval  Research  Laboratory 

ATTN:  Code  6707,  K.  Whitney 
ATTN:  Code  7550,  J.  Davis 
ATTN:  Code  6701 

Naval  Surface  Weapons  Center 
ATTN:  Code  F31 

Strategic  Systems  Project  Office 
ATTN:  NSP 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Geophysics  Lab 
ATTN:  PH,  C.  Pike 

Air  Force  Weapons  Lab 
Air  Force  Systems  Command 
ATTN:  SUL 
ATTN:  NT 
ATTN:  NXS 
2  cy  ATTN:  DYC 

Ballistic  Missile  Office 
Air  Force  Systems  Command 
ATTN:  MNNH 
ATTN:  MNRTE 
ATTN :  KNNG 

Deputy  Chief  of  Staff 
Research,  Development  &  Acq 
Department  of  the  Air  Force 
ATTN:  AFRDQI 

Headquarters  Space  Division 
Air  Force  Systems  Command 
ATTN:  SKF 

Rome  Air  Development  Center 
Air  Force  Systems  Command 
ATTN:  ESR,  E;  Burke 

Strategic  Air  Command 
Department  of  the  Air  Force 

ATTN:  NRI-ST1NF0  Library 
ATTN:  XPFS 


OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 

ATTN:  OSWR/STD/HTB,  A.  Padgett 

NASA 

Lewis  Research  Center 
ATTN:  N.  Stevens 
ATTN:  C.  Purvis 
ATTN:  Library 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

Lawrence  Livermore  National  Lab 

ATTN:  Technical  Info  Dept  Library 

Los  Alamos  National  Scientific  Lab 
ATTN:  MS  364 

Sandia  National  Labs 

Livermore  Laboratory 
ATTN:  T.  Dell  in 

Sandia  National  Laboratories 
ATTN:  3141 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerospace  Corp 

ATTN:  Library 
ATTN:  V.  Josephson 

AVCO  Research  &  Systems  Group 
ATTN:  Library  A830 

Boeing  Company 

ATTN:  P.  Geren 

Computer  Sciences  Corp 
ATTN:  A.  Schiff 

Dikewood  Corporation 

ATTN:  Technical  Library 

Dikewood  Corporation 

Los  Angeles  Operations 
ATTN:  K.  Lee 

EGSG  Hash  Analytical  Services  Ctr,  Inc 
ATTN:  Library 

Ford  Aerospace  &  Communications  Corp 
ATTN:  A.  Lewis 
ATTN:  Technical  Library 

General  Electric  Company 

Space  Division 

ATTN:  J.  Peden 

General  Electric  Company— TEMPO 
ATTN:  W.  McNamara 
ATTN:  DAS 1 AC 

Hughes  Aircraft  Company 

ATTN:  Technical  Library 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued 

Institute  for  Defense  Analyses 
ATTN:  Classified  Library 

IRT  Corp_ 

ATTN:  B.  Williams 
ATTN:  Library 
A1TN:  N.  Rudie 

JAYCOR 

ATTN:  E.  Wenaas 
ATTN:  Library 

JAYCOR 

ATTN:  R.  Sullivan 

Johns  Hopkins  University 
Applied  Physics  Lab 

ATTN:  P.  Partridge 

Kaman  Sciences  Corp 
ATTN:  U.  Rich 
ATTN:  Library 
ATTN:  N.  Beauchamp 
ATTN:  D.  Osborn 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  Dept  85-85 

McDonald  Douglas  Corp 

ATTN:  S.  Schneider 

Mission  Research  Corp 

ATTN:  C.  Longmire 
ATTN:  R.  Stettnner 

Mission  Research  Corp 
ATTN:  B.  Goplen 

Mission  Research  Corp— San  Diego 
ATTN:  V.  Van  Lint 
ATTN:  Library 

Pacific-Sierra  Research  Corp 
ATTN:  H.  Brode 

R  &  D  Associates 

ATTN:  S.  Siegel 
ATTN:  Technical  Info  Center 
ATTN:  L.  Schlessinger 
ATTN:  R.  Schaefer 
ATTN:  P.  Haas 

Rockwell  International  Corp 
ATTN:  Library 

Science  Applications,  Inc 
ATTN:  W.  Chadsey 

Spire  Corp 

ATTN:  R.  Little 

SRI  International 

ATTN:  Library 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continue 


Systems,,  Science  A  Software,  Inc 
ATTN:  A.  Wilson 
ATTN:  Library 
ATTN:  M.  Kandell 
ATTN:  Di  Parks 
ATTN:  G.  Schnuelle 
ATTN:  I.  Katz 
ATTN:  E.  Waisman 
ATTN:  P.  Steen 
ATTN:  M.  Rotenberg 


DEPARTMENT  OP  DEFENSE  CONTRACTORS  (Continued) 

TRW  Defense  S  Space  Sys  Group 
ATTN:  D.  Clement 

ATTN:  Technical  Information  Center 

Hughes  Aircraft  Company 
ATTN:  E.  Smith 
ATTN:  W.  Scott 
ATTN:  A.  Narevsky 


